Using HI to probe large scale structures at z ~ 3 . 
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ABSTRACT 

The redshifted 1420 MHz emission from the HI in unresolved damped Lyman-a 
clouds at high z will appear as a background radiation in low frequency radio obser- 
vations. This holds the possibility of a new tool for studying the universe at high-z, 
using the mean brightness temperature to probe the HI content and its fluctuations 
to probe the power spectrum. Existing estimates of the HI density at z ~ 3 imply a 
mean brightness temperature of lmK at 320 MHz. The cross-correlation between the 
temperature fluctuations across different frequencies and sight lines is predicted to vary 
from 10 -7 K 2 to 10 -8 K 2 over intervals corresponding to spatial scales from 10 Mpc to 
40Mpc for some of the currently favoured cosmological models. Comparing this with 
the expected sensitivity of the GMRT, we find that this can be detected with ~ lOhrs 
of integration, provided we can distinguish it from the galactic and extragalactic fore- 
grounds which will swamp this signal. We discuss a strategy based on the very distinct 
spectral properties of the foregrounds as against the HI emission, possibly allowing the 
removal of the foregrounds from the observed maps. 

Subject headings: cosmology:theory, observations, large scale structures - diffuse radia- 
tion. 
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1. Introduction 

The problem of determining the distribution of matter on large scales in the universe and un- 
derstanding the large scale structure (LSS) formation is of prime importance in modern cosmology. 
Observing the angular positions and redshifts of galaxies has been the most straightforward method 
of mapping the LSS in the present day universe (Peebles 1993, Peebles 1980), but the interpretation 
of these observations is complicated by the fact that the relation between the distribution of galax- 
ies and the distribution of underlying matter is not fully understood (Bardeen et al. 1986). Other 
direct methods use galaxy clusters or super-clusters as tracers to map the large scale structures. A 
large variety of techniques have been developed and applied to quantify the distribution of galaxies, 
and amongst these the correlation functions (two-, three-point correlation functions, etc) and their 
Fourier counterparts (power spectrum, bi-spectrum, etc.) have been most popular (Peebles 1980). 
Much of the work comparing the observations with different theories has been based on these. 

An alternative approach is to use observations of fluctuations in the diffuse extragalactic back- 
ground radiation at different wavelengths to probe the large scale structure. Here, the observations 
of anisotropies in the cosmic microwave background radiation (CMBR) have been most impor- 
tant. These observations probe the clustering of matter at the last scattering surface (e.g. Bond 
1996), and combined with the information from the galaxy surveys, they have been successful in 
narrowing down the allowed class of theoretical models of LSS formation (e.g., Peacock 1999 and 
reference therein) . The study of diffuse background at other wavelengths has been used to infer the 
clustering properties of matter at more recent epochs. Shectman (1974) observed fluctuations in 
the optical background; the results of this observation are consistent with predictions from galaxy 
counts (Peebles 1980). This method has also been applied to far- infrared background (Kashlinsky 
et al. 1997) and recently the first detection of fluctuations in this background has been reported 
(Kashlinsky et al. 1999). Gorjian, Wright and Chary (1999) have recently reported a tentative 
detection of a near-infrared background. There are similar predictions for fluctuations in the X-ray 
background (Barcons et al. 2000, Barcons et al. 1998). 

In this paper we investigate the possibility of using the extragalactic background radiation at 
low radio frequencies (meter wave) to probe the LSS. This is largely motivated by the fact that the 
Giant Meter-wave Radio Telescope (GMRT; Swarup et al. 1991) which is designed to observe in 
several frequency bands in the interval 150 MHz to 1420 MHz has recently started functioning. This 
frequency range corresponds to the 1420MHz emission from HI in the redshift range < z < 8.5. 
Absorption studies along sight lines to quasars indicate that most of the HI in the redshift range 
< z < 3 is in damped Lyman-a (DLA) clouds and the density of HI in DLAs has been determined 
reasonably well from absorption studies (Lanzetta, Wolfe, & Turnshek 1995). Although the HI 
emission from individual DLAs at high z is too faint to detect using presently available telescopes, 
the redshifted HI emission from unresolved DLAs will contribute to the background radiation at low 
frequencies. In this paper we investigate how the brightness temperature of this radiation is related 
to the density and peculiar velocity of the HI. We consider the possibility of detecting correlations 
in the fluctuations in this component of the background radiation and using this to probe LSS at 
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high z. 

The possibility of observing the redshifted 21 cm emission from HI at high z has been discussed 
earlier by many authors in a large variety of contexts. One of the first papers (Sunyaev & Zel'dovich, 
1975) considers the possibility of meter-wave observations to detect protogalaxies and protoclusters 
at z < 10. There have been several attempts to detect the HI in proto-clusters and proto-super 
clusters (see Subrahmanyan & Anantharamaiah 1990 and reference therein). 

More recently Subramanian and Padmanabhan (1993) have calculated the abundance of pro- 
toclusters as a function of the redshifted HI flux density for various redshifts for both the CDM and 
HDM models. In a later paper Kumar, Padmanabhan and Subramanian (1995) have calculated 
the line profile of the HI emission from a spherically symmetric protocluster. Bagla, Nath and Pad- 
manabhan (1997) and Bagla (1998) have used a combination of N-body simulations and a model 
for the behaviour of the baryons to calculate the abundance and the expected flux from the HI in 
structures like protoclusters at high redshifts. A major uncertainty in all of these works is in the 
assumptions about the HI content of the universe at high redshifts. The main focus of all of these 
works has been on individual peaks of the density fluctuations (protoclusters) which will manifest 
themselves as detectable features in low-frequency radio maps. Subramanian and Padmanabhan 
(1993) have also studied the possibility of detecting the excess variance in radio images due to the 
density fluctuations in the HI at high z. 

Katz, Weinberg and Hernquist (1998) have used smoothed particle hydrodynamic simulations 
to study the distribution of HI gas at high z and they consider the possibility of detecting HI 
emission from galaxies at z > 2. 

Tozzi et al. (1999) and Shaver et al. (1999) have studied the possibility of detecting the HI 
in the IGM at z > 5. The state of the HI at these redshifts is unknown and these works are based 
on different scenarios for the reionization of the universe. 

The work presented here differs significantly from the previous papers in that: 

(1) It is restricted to z < 3.5 where the HI content of the universe is well determined from DLA 
absorption studies. We use the results of these observations as inputs to our calculations. 

(2) Instead of looking at the possibility of detecting individual features (as has been the focus of a 
large number of previous papers) we have studied the statistical properties of the fluctuations 
in the brightness temperature in low frequency radio maps. The statistical quantity we 
have studied is the cross-correlation between the temperature fluctuations along different 
lines of sight in radio map made at different frequencies. Individual features corresponding 
to protoclusters are rare events and protoclusters with flux in the range 1.5 — 3mJy are 
predicted to occur with abundances in the range 10" 8 - 10" 7 Mpc~ 3 in the CDM model 
(Subramanian and Padmanabhan, 1993). Even small density fluctuations will contribute to 
the cross-correlation signal and our proposal has the advantage of simultaneously using the 
signal in all the pixels in all the frequency channels across the bandwidth of observation. 



-4 - 



The structure of the paper is as follows: in §2 we discuss the basic formalism of calculating 
the brightness temperature and fluctuations in the redshifted 21 cm radiation, in §3 we discuss 
the numerical results for two currently favoured cosmological models. In §4, the observational 
difficulties in the presence galactic and extragalactic foregrounds are presented, and we discuss a 
possible strategy for overcoming these. §5 gives a summary of our main results. 



2. Formalism. 

We treat the HI in damped Lyman-a (DLA) systems as a continuous distribution with nm(x, t) 
denoting the comoving number density of HI atoms in the excited state of the hyperfine transition. 
Such a treatment is justified in a situation where the resolution of the observations is not sufficient 
to detect individual DLAs. In addition, the fact that the HI actually does not have a continuous 
distribution but is distributed in discrete objects can, if required, be taken into account when 
calculating statistical measures of the fluctuation in the HI distribution. 

The HI emission which is at a frequency v e = 1420MHz in the rest frame of the gas it is 
emitted from is redshifted to a frequency v for an observer located at the origin of the coordinate 
system. Taking into account the effects of both the expansion of the universe as well as the peculiar 
velocity v(x,t) of the HI at the time of emission and t>(0,io) of the observer at the present time, 
the comoving coordinate x of the HI and v are related by 



x = nc I 2 1° \ (1) 

J t v t ^ a z H(a) 



•Ai-w) 



in a spatially flat universe. Here n = x/x is a unit vector along the line of sight to the HI, 
W(x) = n ■ [v(x,t) — v(0,t(j)]/c accounts for the effects of the peculiar velocities and 

a(t) 



H(a) = = HoVn^a-z + A (2) 

is the Hubble parameter at the epoch when the scale factor has value a. 

When discussing observations of the redshifted HI emission, it is convenient to use v = nu to 
simultaneously denote the frequency and the direction of the observation. The vector v fixes both 
the comoving position x of the HI from where the radiation originates and the time t at which the 
radiation originates, and we shall use v and (x,t) interchangeably. Here we calculate how T(v) 
the brightness temperature of the radiation is related to the density nE_i(x,t) and the the peculiar 
velocity v(x,t) of the HI. 

The energy flux in the frequency interval d 3 v can be calculated if we know the number of 
excited HI atoms in the comoving volume d 3 x from where the radiation originates and it is given 
by 

h P i/ e A 2 in m (x,t)d 3 x 
Flux = T — (3) 
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Here tl is the luminosity distance which is given by 



r L (x) = x[l-W]-^ 



The flux is also related to the specific intensity as follows 



(4) 



Flux = 



I{u) 



d 6 u 



which allows us to calculate the specific intensity to be 

I{v) _ h P is e A 21 n H i(u) 



1 


f 


dx 


W{u)] Xj 




du 



(5) 



(6) 



where is the Jacobian of the transformation from u to x given in equation (1). 

It should be noted here that the effect of the peculiar velocity (i.e.W) can be neglected if we 
restrict the analysis to scales which are much smaller than the horizon. This is not true for the 
spatial derivatives of the peculiar velocities which appears in the Jacobian and we retain such terms 
in our analysis. Calculating the Jacobian gives us 



dx 



du 



/ x \ 2 c v, 



1 



v e (n • V)(n • v) 
v H(^) 



(7) 



Equations (6) and (7) allows us to calculate the brightness temperature of the radiation 

" 2 N /(„) 



1k B 



which gives us 



T(u) = 



T{u)- 
T 21 N 2 i(u) A 21 (v, , ■ 



8vr 



H( 



^3, 



1 



v e (n ■ V)(n • v) 



(8) 



(9) 



where T 2 \ = h p u e /kB and N 2 \(u) = {c/v e ) nui{u) is the number of HI atoms in the excited state 
in a comoving volume (21cm) 3 . 

The number density nu\{u) can be written as a sum of two parts, namely the mean iIhi(u) 
and the fluctuation Are#/(i/). We assume that the fluctuation in the number density of HI can 
be related to the perturbations in the underlying dark matter distribution S(u) through a time 
dependent linear bias parameter b(u) which gives us 



N 21 (u)=N 21 (u)[l + b{u)5{u)}. 

We use this to calculate the isotropic part of the temperature 

2.38 K/t- 1 jV 2 iH /z/ e \ 2 



T{y) 



(10) 



(11) 
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and the fluctuation 



AT(iz) = f{u) 



b{v)5{v) - 



u e (n • V)(n • v) 



(12) 



Now, in the linear theory of density perturbations, both the perturbation and the peculiar 
velocity can be expressed in terms of a potential i.e. 



S(x,t) = D(t)V 2 ip(x) 



(13) 



and 



v(x, t) = -f{VL m )H{t)a{t)D{t)V^{x) . (14) 

where D{t) is the growing mode of linear density perturbations (Peebles 1980). In a spatially flat 
universe the function / can be well approximated by the form (Lahav et. al. 1991) 



/(Or 



where the time dependence of Q m can be expressed as 



Vl m (v) 



H(-) 



n 



m0 



(15) 



(16) 



Using this and defining a time dependent linear redshift space distortion parameter (5{y) = f(Q m ) /b{y) 
we can express the fluctuation in the temperature as 



A2» = T A {y)b{y)D{y) [V 2 + (3(u)(n ■ V) 2 ] ^{x) (17) 

We next calculate the cross-correlation between the temperature fluctuations along two differ- 
ent lines of sight ni and ri2 at two different frequencies v\ and z^. The quantity we consider is the 
correlation 

w{v u u 2 ) = (AT(i/i)AT(i/ 2 )> (18) 

which is a function of the two frequencies v±, v<i and 9 the angle between the two lines of sight. 
Using equation (17) and defining (j)(xi2) = (^(^1)^(^2)) to be the two point correlation of the 
potential ip(x), we obtain 

w{u u v 2 ) = f x f 2 D x D 2 b x b 2 [V 2 + ^(ni • V) 2 ] [V 2 + (3 2 (n 2 ■ V) 2 ] <f>(x 12 ) (19) 

where x\ 2 =| X\ — x 2 \ and we have used the notation T\ = T{v{), etc. 

The function <p(x) is related to the correlation of the perturbations in the underlying dark 
matter distribution, and ti, x 2 , t 2 ) = (8(xi, t\)5(x 2 , t 2 )) the two point correlation of the per- 
turbation in the dark matter density at the point xi at the epoch t\ and x 2 at the epoch t 2 can be 
written in terms of the potential <fi as 

h, x 2 , t 2 ) = D(h)D(t 2 )i(x 12 ) = £>(ti) J D(t 2 )VV(^i2) • (20) 
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We have introduced the function £(£12) so that we can write £(£Ci, t±, X2, £2) as a product of two 
parts, one which has the temporal variation and another the spatial variation. Equation (20) can 



be inverted to express the different derivatives of 
moments of £(#12) which are defined as 



(X12) which appear in equation (19) in terms of 



L(X12) 



n + 1 T 12 z 



r n+l 
J 12 



t(y)y n dy. 



(21) 



The form of the angular correlation is further simplified if we restrict 9 to be very small. Under 
this assumption n\ ~ ri2, and we use n to denote the common line of sight. We also use \i = 
n ■ {x\ — X2)jx\2 for the cosine of the angle between the line of sight n and the vector x\ — X2 
joining the two points between which we are measuring the correlation. The relation between the 
different distances and angles is shown in figure 1, and we have x\2 = y/x'f + x\ — 2x\X2 cos(#) and 
H = (xi — X2cos{0))/x\2- Using these we finally obtain the following expression for the two point 
correlation of the temperature 



w(v 1 ,v 2 ,9) = T 1 T 2 i}iZW2{[(l + /V)(l + /V)]£>i2) 



+ 



(I - V 2 )(Pi + ft) + {\~ 3m 2 + 



6(^12) 



10 



7 



£4(2:12) 



(22) 



3. Predictions for different models. 

The density of HI (f2m) in DLAs has been determined reasonably well for < z < 3 from 
absorption studies (Lanzetta, Wolfe, & Turnshek 1995) and they find that the observed evolution 
of Qui is well approximated by Clui(z) = OHioexp(az) with Ohio = 0.18 ± 0.04 x 10 _3 /i _1 and 
a = 0.60 ± 0.15 for q = 0, and Ohio = 0.19 ± 0.04 x lO" 3 /^ 1 and a = 0.83 ± 0.15 for q = 0.5, 
We have used this to calculate the isotropic part of the background brightness temperature due to 
the redshifted HI emission (equation 17). We have considered two spatially flat FRW cosmological 
models with parameters (I) f2 m o = 1.0 and Ao = 0, and (II) VL m Q = 0.3 and Ao = 0.7, with h = 0.5 
in both cases. It should be pointed out that the fits given Lanzetta et al. (1995) are not valid for 
the model with a cosmological constant, and we use the qo = fit in this case. 

Figure 2 shows the brightness temperature as a function of frequency for both the cases. 

We find that the temperature increases rapidly as we go to lower frequencies (higher z) and 
it is around 1 mK at v ~ 330MHz which corresponds to z ~ 3. The increase in the temperature 
is a direct consequence of the increase of HI density with increasing redshift. The HI density is 
not very well determined at z > 3 and there is evidence that Qm falls off at higher redshifts 
(Storrie-Lombardi et al. 1996). 

We next consider w(i/i,i>2,9) the cross-correlation in the temperature fluctuation at different 
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frequencies which we have calculated for the two cases discussed above. We have calculated the 
dark-matter two point correlation function (equation 20) using the analytic fitting form for the 
CDM power spectrum given by Efstathiou, Bond and White (1992). For case (I) we use a value 
of the shape parameter T = 0.25 and for (II) we use T = 0.14. The power spectra are normalized 
using the results of Bunn and White (1996) based on the 4-year COBE data. It should be noted 
that he redshift evolution of the matter correlation function and f(fl m ) are quite different in the 
two cases that we have considered. The predictions for w(vi,i>2-,9) are shown in figures 3 and 4. 
We have assumed that the HI faithfully traces the matter distribution and set b(u) = 1 throughout. 

In our analysis we have kept v\ fixed at 320MHz and let vary over a band of 16MHz 
centered around v\, while 9 takes values upto 2°. In order to measure w(ui,V2,9) over this range 
we would need radio images of the background temperature fluctuation in a 2° x 2° field at different 
frequencies in a 16MHz band centered around 320MHz and w(ui,V2,9) could be estimated using 

w{u l ,u 2 ,9) = (AT U1 (ni)AT U2 (n 2 )) (23) 

where AT ul (ni) refers to the temperature fluctuation along the direction n\ in the image made at 
frequency v\ , and the angular brackets denote average over all pairs of directions n\ and n 2 which 
are separated by an angle 9. The central frequency, bandwidth and angular range have been chosen 
keeping in mind the Giant Meter wave Radio Telescope (GMRT, Swarup et al. 1991) which has 
recently become operational. The frequency intervals and angular separation can be converted to 
a corresponding comoving length-scale and in case (I) 1MHz = 8.9Mpc and l' = 1.8Mpc, while in 
case (II) 1MHz = 16.0Mpc and l' = 2.7Mpc. The linear theory of density perturbations which we 
have used here can be applied at scales > lOMpc, which covers most of the region shown in the 
figures. 

We find that w{yi,U2,9) is between 10 _7 K 2 to 10 _8 K 2 when the comoving-distance corre- 
sponding to the separation in frequencies and direction is less between lOMpc and 40Mpc beyond 
which it falls off. The cross-correlation between the temperature fluctuations is expected to be 
larger at small scales where the linear theory of density perturbations will not be valid and we 
have not considered these scales here. The temperature fluctuations are anti-correlated when the 
separation in frequencies exceeds the angular separation. This occurs because of the effect of the 
peculiar velocity which produces a "distortion" very similar to the effect it has on the two-point 
correlation function in redshift surveys. 

4. Observational prospects 

In this section we discuss the prospects of actually observing the background radiation from 
the HI in damped Lyman-a clouds. Our discussion is restricted to observations at around 320 MHz 
largely because we have reliable estimates of the HI content at z < 3 and there are indications that 
z ~ 3 might be the redshift where the HI content of the universe is maximum (Storrie-Lombardi et 
al. 1996). Another reason for our choice of this frequency range is that the GMRT (Swarup et al. 



-9- 



1991) which is already functioning at these frequencies is expected to have an angular resolution 
of 10" and reach noise levels of 100/t/K in around 10 hrs of integration. The angular resolution 
and sensitivity of the GMRT will be sufficient for detecting both the isotropic component as well 
as the correlations in the component of the background radiation arising from the HI at z — 3 
(Figures 2, 3 and 4) provided we can distinguish this component from the contribution due to other 
sources. 



4.1. Galactic and extra-galactic foregrounds 

Any observations at low frequencies will have a very large contribution from the synchrotron 
radiation from our own galaxy. Observations at 408 MHz (Haslam et al. 1982) with resolution 
l°xl° indicate a minimum temperature of ~ 10 K at 408 MHz. Using T oc z/~ 2 - 7 as indicated by 
the spectral index of galactic synchrotron radiation, the temperature at v ~ 320 MHz is ~ 20 K. 
Comparing this number with the expected background from redshifted HI radiation (Figure 2) the 
galactic foreground is seen to be several orders above the expected signal. Another quantity of 
interest to us is the fluctuations in the galactic radiation at angles ^,1°. However as the resolution 
of Haslam maps is ~ 1° it cannot be used to make any predictions about the fluctuations at the 
angular scales of interest. 

Another source of contamination is the continuum radiation from unresolved extragalactic 
sources (the resolved ones will be removed from the image before analysis) and we use results from 
the recent FIRST survey (for details see White et al. 1997 and references therein) to estimate this. 
Tegmark & Estathiou (1995) provide an analytic fit for the number of sources per unit flux 0[mJy] 
per steradian in this survey, and at flux levels <j> <C lOOmJy this can be approximated by 

dn _ 5.24 x 10 5 / cf> \ _1 ' 65 
dcf) mJysr \0.75mJyy 

We use this to estimate the total contribution from sources fainter than 100 /itJy, and converting 
this to brightness temperature we find that the continuum emission from unresolved radio sources 
is expected to produce a background radiation with brightness temperature ~ .IK at 1.5 Ghz. For 
the purposes of estimating an order-of-magnitude we assume this value to be representative of what 
we expect at ~ 320 MHz. 

The two-point correlation function of the sources detected by FIRST has been estimated 
(Cress et al. 1996): w(9) ~ O.180 -11 , with 9 in arcminutes. Little is known about the correlations 
expected in the fluctuations in the contribution to the background radiation from unresolved radio 
sources. 

It is clear from the foregoing discussion that both galactic and extragalactic continuum signals 
are likely to be so high that they would totally swamp the signal which we want to detect and 
therefore it is not possible to directly detect the radiation from HI unless we find some method for 
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removing the foregrounds. 

At present the most promising strategy is to use the proposed observations themselves to 
determine the foreground and remove this from the data, and we next discuss a possible method 
for doing this. 



4.2. Removing foregrounds 

We consider the GMRT band centered at 325 MHz for our discussion in this section. The total 
bandwidth of 16 MHz is divided in 128 channels each with Ai> = 1.25 kHz, and the signal in any 
frequency channel is expected to be dominated by the foregrounds. 

Here we briefly review the standard method of continuum subtraction which is generally used 
in spectral line detections and then discuss how this method can be used for the analysis proposed 
here. The reader is referred to Subrahmanyam and Anantharamaiah (1990) and references therein 
for examples of how this method is applied in searches for HI emission from protoclusters. 

We represent the observed signal yi(n) in the i th frequency channel in the direction n as: 

Vi {n) = Xi {n) + fi{n) + Ni (25) 

where Xj(n) and fi(n) are the contribution to the output signal from the HI and the foregrounds 
respectively, and Ni is the receiver noise. 

We should at this stage remind ourselves that foregrounds are sources which are emitting 
continuum radiation — the primary such source being the synchrotron radiation from our own galaxy. 
We also have the contribution from unresolved extragalactic radio sources which, as we have seen 
above, is a smaller contribution. The sum total of these along any line of sight is expected to be 
a smooth function of the frequency. On the contrary the contribution from the HI will come from 
individual damped Lyman a clouds along the line of sight with a velocity width ~ 200 km sec -1 
(Prochaska & Wolfe 1998) and this will correspond to lines of width: 

Ay 

Au ~ — v ~ 0.2MHz. (26) 

Although the total number of damped Lyman-a clouds in a 1° x 1° x 16 MHz field is expected 
to be quite large (~ 3 x 10 4 ) we expect to have only one damped Lyman-a cloud along a single 
20 x 20 x 16 MHz synthesized GMRT beam. The probability of two damped Lyman-a clouds 
occurring right next to each other in the same synthesized beam and causing confusion is quite 
small, and hence the signal we are looking for will come in the form of lines of width ~ 0.2MHz 
which we will not be able to detect as the individual lines will be swamped by both the foreground 
as well as the noise. What we can do is to fit the smooth component of the signal yi(n) by a 
function Fj(n) and subtract this from the output signal to remove the effect of the foregrounds 

Si(n) = yi (n) - Fi(n) = x;(n) + A^n) + E^n) (27) 
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where we call Si(n) the reduced signal and Ei(n) is the possible error in the fitting procedure. If 
our foreground subtraction works correctly, the reduced signal should have contributions from only 
the HI signal and the noise. 

We first determine the mean reduced signal 

Si = (Si(n)) = (Xi(n)) + {Ni{n)) + (E^n)} (28) 

where the angular brackets denote average over all lines of sight. If the errors in the foreground 
subtraction can be made smaller than the signal i.e. {Ei(n)) < (xj(n)) then this will give an 
estimate of the background brightness temperature due to the HI emission i.e. T{y) = Si where 
the index i refers to the channel with central frequency v. 

We next consider the quantity 

ASi(n) = Si(n) - S t (29) 

which gives the fluctuation in the brightness temperature. We can use this to estimate the cross- 
correlation in the fluctuations at different frequencies 

w{yx, vz, 9) = (AS^ASjim)} (30) 

where i and j are the channel with frequency v\ and respectively, and n and m are two lines of 
sight separated by an angle 9 and the angular brackets denote average over all such pairs of lines 
of sight. 



5. Summary and Conclusions 

We have investigated the contribution from the HI in unresolved damped Lyman-a clouds 
at high redshifts to the background radiation at low frequency radio waves (meter waves). The 
isotropic part of this radiation depends on the density of HI and the background cosmological model, 
while the fluctuations in this component of the background radiation have an added dependence 
on the fluctuations in the distribution of the damped Lyman-a clouds and their peculiar velocities. 

We have used estimates of the HI density available from absorption studies to calculate the 
brightness temperature of this radiation. We find that this has a value ~ 1 mK at 320 MHz which 
corresponds to z ~ 3. The distribution of damped Lyman-a clouds is assumed to trace the un- 
derlying dark matter distribution which also determines the peculiar velocities. Using this and the 
linear theory of density perturbations, we have calculated the relation between the fluctuation in 
this component of the background radiation and the density perturbations at high z. Observations 
of the cross-correlations of the fluctuations at different sight lines across images produced at dif- 
ferent frequencies holds the possibility of allowing us to probe the two point correlation function 
(or power spectrum) at high redshifts. We have calculated the expected cross-correlations for two 
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currently acceptable CDM models and find it to be in the range 10 _7 K 2 to 1CT 8 K 2 at v ~ 320MHz 
for separations in sight lines and frequencies such that the corresponding spatial separation is in 
the range lOMpc and 40Mpc. The cross-correlations are expected to be larger at smaller scales 
where the linear theory cannot be applied. Our results show that both the isotropic background 
(Figure 2) and its fluctuations (Figure 3 and 4) can be detected by GMRT which is the largest 
telescope operating at meter waves at present, provided this signal can be distinguished from other 
sources which contribute to the low frequency background radiation. 

The biggest obstacles in detecting the HI contribution are the galactic and extra-galactic 
foregrounds, both of which are many orders larger than the signal we want to detect. The fact 
that both those sources of contamination emit continuum radiation while the HI contribution is 
from individual damped Lyman-a clouds each of which emits a spectral line with a relatively 
small velocity width keeps alive the possibility of being able to distinguish this signal from the 
contamination. We have, in this paper, considered one possible approach which might allow us to 
model and subtract the foreground along any line of sight. More work is needed in this direction 
and work is currently underway in investigating other viable possibilities for foreground removal. 

All the authors would like to thank Jayaram Chengalur for many useful discussions on issues 
related to the GMRT and foregrounds. 
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Fig. 1. — This shows the various distances and angles involved in calculating w(y\, ^2,0), where v\ 
can be converted to x\ using equation (1). 
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Fig. 2. — This shows the isotropic part of the brightness temperature at different frequencies. The 
solid line shows model (I) and the dashed line model (II) of section 3. 
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Fig. 3. — For Q m o = l,h = 0.5, (a) shows w(v\,V2-,0) in K 2 vs. 6 with v\ = 320MHz and curves 
(1),(2) and (3) corresponding to v 2 = 320MHz, 322MHz and 324MHz respectively; (b) shows 
contours of equal w(u\, f2,0) at logarithmic intervals of w with v\ = 320MHz. Here 1MHz = 8.9Mpc 
and l' = 1.8Mpc. The dashed lines show negative values of w. 
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Fig. 4.— For Q m0 = 0.3, An = 0.7, h = 0.5, (a) shows w{v 1 ,v 2 ,0) in K 2 vs. 9 with v x = 320MHz 
and curves (1),(2) and (3) corresponding to v 2 = 320MHz, 322MHz and 324MHz respectively; 
(b) shows contours of equal w(i/i,i>2,8) at logarithmic intervals of w with v\ = 320MHz. Here 
1MHz = 16.0Mpc and l' = 2.7Mpc. The dashed lines show negative values of w. 



